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Abstract 

We investigate possible parafermionic states in rapidly rotating ultracold bosonic atomic gases at 
lowest Landau level filling factor v = fc/2. We study how the system size and interactions act upon 
the overlap between the true ground state and a candidate Read-Rezayi state. We also consider 
the quasihole states which are expected to display non-Abelian statistics. We numerically evaluate 
the degeneracy of these states and show agreement with a formula given by E. Ardonne. We 
compute the overlaps between low-lying exact eigenstates and quasihole candidate wavefunctions. 
We discuss the validity of the parafermion description as a function of the filling factor. 

PACS numbers: 03.75.Lm, 03.75.Kk, 73.43. Cd, 73.43.Nq 
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I. INTRODUCTION 



Rotating Bose-Einstein condensates display a wealth of interesting physics. One of the 
most striking achievements in this field is the observation of the Abrikosov lattice of vor- 
tices^. With increasing rotation speed, it has been predicted that this lattice will melt and 
is replaced by more exotic quantum phases. When the rotation frequency is close to the har- 
monic trapping frequency and strong confinement is applied along the rotation axis, strongly 
correlated states belonging to the family of quantum Hall liquid states should appea r 3 i 4 i 5 i 6 . 

Here we study bosonic atoms with only one hyperfine species (i.e. spinless bosons) in 
such regime. We assume that the temperature is low and the interactions are weak enough 
so that the lowest Landau level (LLL) approximation is valid. The system may then display 
the fractional quantum Hall effect (FQHE) as in two dimensional electron systems (2DES) 
under strong magnetic field. The Coulomb interaction is replaced by the s-wave scattering 
between the ultracold atoms. An analog of the filling factor v for 2DES can be defined : 
indeed v = N/N^ is the ratio between the number of atoms N and the number of vortices 

that would be present in the system if it was a Bose condensate. The quantity 
is the equivalent of the number of flux quanta in 2DES systems. In this regime, several 
fractions have been predicted. The most prominent one 3 appears at u = 1/2, for which 
the Laughlin state is the exact ground state 7 . Evidences for other fractions from the Jain 
principal sequence v = p/(p + 1) such as v = 2/3 and v = 3/4 have been pointed out&& and 
can be understood within the composite fermion theory^. 

Hierarchical quantum Hall states are not the only interesting states that have been pre- 
dicted below the critical filling factor where the lattice of vortices mell£. Due to the bosonic 
statistics, if we assume that the equivalent cyclotron gap is large enough, we can have filling 
factors greater than one and still stay entirely in the LLL. Within this hypothesis, even 
more exotic states should appear for fractions v = k/2. The first one is the Moore- Read 
(MR) stated (or Pfaffian state) that should occui^ at v — 1 (k = 2). This state was first 
introduced to explain the fermionic fraction v = 5/2 in 2DES. Higher k values are associ- 
ated to the so-called Read-Rezayi (RR) states*^. Because of the parafermionic behavior 
of these states, their excitations have surprising non-Abelian statistics. So far, there are no 
well-established physical situation where these states play a role. The original suggestion by 
Read and Rezayi is that they may explain the incompressible states observed in the second 
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Landau level on the flanks of the elusive v = 5/2 state. 

The RR states (or clustered states) in ultracold rotating atomic gases have been already 
the focus of several works. Numerically exact diagonalizations of small systems have provided 
some hints of the presence of RR states. In the seminal work by Cooper, Wilkin and 
Gunn^, spectra in the torus geometry exhibit the special ground state degeneracy associated 
with the topology of the RR states and have excellent overlaps with the explicit RR trial 
wavef unctions. In the spherical geometry there is also a set of incompressible states with 
the special relationship between the flux and the number of particles of the RR states 6 . 
Extrapolation of the gap points to a non-zero value for the MR v = 1 case, whereas the 
v = 3/2 and v = 2 results do not show clear evidence for a smooth thermodynamic limit. 
On the sphere geometry, the overlap is excellent for the MR state and tend to a nonzero 
value as the system size increases^. A more recent work 12 has been done on the v = 3/2 
case. It shows that the overlap can be improved by adding a longer range dipole-dipole 
interaction 

Our purpose is to go beyond existing studies and look at size effects for fractions 
v = 3/2,2,5/2 using exact diagonalizations on the sphere. We also check if the quasihole 
ground states are present at these filling factors by evaluating overlap between subspaces 
spanned by these states and the lowest energy excitations of the delta function interaction, 
s-wave scattering system. Appearance of such quantum states with the correct degeneracy 
predicted by conformal field theory arguments is a strong hint of the validity of the RR state 
hypothesis. 

In section [HJ we give an overview of the clustered states and their excitations while 
section ITTT1 is devoted to the conformal field theory (CFT) formulation. Section HVl is a brief 
description of the numerical method we use. In section we give the results for the overlap 
of the ground states. We discuss how the system size, longer range or higher order n-body 
interaction impact on the overlap. Section IVII is devoted to the quasihole excitations. In 
addition to the overlap values, we also give numerical evaluation of quasihole degeneracy on 
sphere for fractions v = 1,3/2, 2, 5/2 and compare them to a formula due to Ardonne^ to 
check the validity of the conformal field theory approach. 
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II. PARAFERMIONIC STATES 



For the sake of simplicity, we use the disk geometry in this section. In the symmetric 
gauge, the LLL one-body wave functions are given by : 



' ,z m e-^ 2 '\ (1) 



Ttlb\ r J t 

v / 27r2 r "m! 

where z = x + iy and we take the magnetic length Ib to be equal to unity. Any iV-body 
wave function of particles in the LLL can be written as a polynomial V in the particle zi 
coordinates : 



^( Zl ,...,z N ) = V(z 1 ,...,z N ) e-^\ 2 " 12 . (2) 

From now on, we drop the global Gaussian factor. The fc-type RR state is the exact zero 
energy ground state of the pure (k + l)-body 5-function interaction hamiltonian : 



Wf B = £ S^( Zil -z i2 )...S^(z ik -z ik+1 ). (3) 

h<...<i k+1 

The corresponding wave function can be written^ : 

^k R = E II X{z<r{kr+1), -,Z a (kr+k)] 
a 0<r<s<N/k 

(4) 

where 

X \ Z X: •■■) z k] Zk+1, z 2k) = {z\ — Zfc+l) 
\Z\ — z k+2) ( z 2 ~ z k+2) ( z 2 — z k+z) 

••• {Zk - z 2k) i z k ~ z k+l) ■ (5) 

The sum is over all permutations a of N elements such that cr(l) < u{k) < ... < a(N — k + 1). 
The number of particles N must be a multiple of k. 

These states are also referred to as clustered states because the wavefunction (jlj) vanishes 
when k + 1 or more particles are at the same position. The /c-type RR state is associated to 
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the filling factor v = k/2. Each RR state is the zero-energy ground state of its corresponding 
Hamiltonian with the smallest total angular momentum. 

The simplest case k = 1 corresponds to the usual Laughlin wave function : 

^Laughlin = J^J ( z i ~ z jf (6) 

and is the exact ground state for rotating bosons with s-wave scattering at v = 1/2 whose 
effective Hamiltonian is given by : 

H L ll = gl 2 B ^2 ^ (2) ( ri ~ r ^ and g = fovcy-- (7) 

where a s is the s-wave scattering length, l z is the characteristic length of the z axis oscillator 
which is used for 2d confinement, and uj c is the cyclotron rotation frequency. 
The case k = 2 is the so called MR/Pfafnan state. It can be rewritten as : 

Vpfaffian = ^ ( ) II ~ ^ ' ( 8 ) 

\Z{ Zj I 

where Pf stands for the pfaffian defined as : 

(1)<t(2) ^a(3)<r(4) • • • At(7V-1)<t(A0 , (9) 

where A is a skew-symmetric N x N matrix (N even), the sum runs over all permutations 
of the index with N values and e a is the signature of the permutation. 

If we deviate from the clustered state at filling factor k/2 by adding A$ vortices (or flux 
quanta in the 2DES analog), quasihole excitations are generated. For each added vortex, k 
quasiholes are nucleated. For the Laughlin state, quasihole ground state wave functions can 
easily be obtained. Any function of the form : 

^Laughlin = P i Z li ~> Z n) * Laughlin, (10) 

where P is a symmetric polynomial, corresponds to a zero-energy many quasihole state. 
For one quasihole at position wi, the polynomial P is just Ylii 2 * ~ w i)- Read an d Rezayi 
have also obtained an explicit formula in the case of the MR stated for two quasiholes at 
positions Wi and w 2 : 

with / (zi, Zj; Wi, w 2 ) = — wi) (zj — w 2 ) + ( z i — w 2 ) (zj — wi). In the general case, the 
quasihole ground states can be written down using the CFT formulation. This formalism 
also reveals their non-Abelian statistics. 
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III. CFT APPROACH 



There is an elegant way to introduce RR states involving CFT—. The key ideaM is to 
express the wave function as a correlator using the algebra of the parafermions^. This 
algebra is defined a set of field {&i(z), ^?k-i} obeying the following operator product 
expansion (OPE) : 

i/>i(z)M*) ~ <W*-^r (A,+A, '~ Al+i,) 

xil> l+l ,{z') + ■ ■ ■ (/ + /'< k), (12) 

-v ( z ~ z ' 

xVw(/) + --- (!'<i), (13) 
MWHz') ~ (z-zT 2Al 

x (l+^(^-/) 2 T(^) + •••), (14) 
(z - z') 

+^j,mz) + ---- as) 

where ip\ = ipk-i, T(z) is the stress-energy tensor, A; is the conformal weight of the field 
ipi, c is the theory central charge and dijt are numerical coefficients. The algebra of 
parafermions corresponds to the choice A; = l(k — l)/k leading to the central charge c = 
2(k — l)/(fc + 2) and to uniquely determined d^y coefficients. 
Read and Rezayi have shown that the following wave function : 

*r cFT = (v>i (zi) --^i m) n c* - h) 2ik , (is) 

i<j 

is equivalent to expression One can easily show that this expression vanishes quadrati- 
cally as k + 1 particles go to the same point using the OPE rules above (|12H15|) . Within this 
formalism, it has been argued that the zero energy quasihole states can be built by inserting 
a spin field for each quasihole into the correlator of (|16j). For k = 2, this spin field o is 
equivalent to the magnetization operator of the Ising model. In the Ising case, the fusion 
rules are given by : 

cr(z)Mz') = ~ l —n-M z ')i ( 17 ) 

[z — z') 1 

a(z)a(z') = 1 I +(z- z'f 8 ^{z>). (18) 
[z — z') ' 



For n = 2A$ quasiholes, the candidate state is then : 

^pfaffian = 00 »4l (z N ) OO -° On)) 

x n^-^) nn^-^) i/2 - ( i9 ) 

i<j i p=l 

The fusion rule (jl8j) leads to a non-trivial degeneracy of the quasihole states : there are 
2™/ 2-1 ways to fuse the spin operators leading to a non zero correlator, thus giving as many 
different wavefunctions. This so-called intrinsic degeneracy is the key of non-Abelian statis- 
tics : exchanging two quasihole coordinates of a given quasihole state will result in a linear 
combination of states of the same family instead of an overall multiplicative phase factor. In 
the case of the spherical geometry that we will discuss later, an additional (extrinsic) degen- 
eracy arise from the Laughlin-like part of Eq. (j20J) . Determining the multiplet decomposition 
of quasihole states in such a case is a challenging tas k 13 ' 15 ' 17 and constitutes a non-trivial 
check of the CFT approach when compared to numerical calculations. More details will be 
given in section fVTl 

For k > 2, the spin field that we have introduced has to be replaced by one of the primary 
field operators of the Zj, parafermion algebra. The guess is to use the operator o\ which 
minimizes the charge of the quasiholes. The wavefunction (j!9|) can then be generalized to : 

^k Qh = (^1 (21) •••'01 {z n )<t x Oi) On)) 

n 

x Uizi-zjf" nn^-^) i/fc - ( 2 °) 

i<j i p=l 

The fusion rules involving <j\ are more comple x 16 ! 18 but the same remarks as for the k = 2 
case apply, meaning they lead to non-Abelian statistics. Notice that for k > 3, such states 
have been proposed to be a robust way to implement quantum computation^. 

IV. NUMERICAL METHOD 

We use exact diagonalizations to study if the RR states are relevant to the physics of the 
fast rotating boson gases at filling factor v = k/2. Numerical calculations can be done on 
various geometry such as the disk, the torus or the sphere. The disk geometry is plagued 
by edge effects and thus closed geometries are preferred when dealing with bulk properties. 
In this paper, all calculations are done on the spherical geometry^. Due to the SU(2) 
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symmetry, states can be classified with respect to their total angular momentum L and its 
projection along one axis L z . Solutions of the one-body problem are given by the monopole 
harmonics^ (a generalization of the spherical harmonics), which take the following form in 
the LLL : 

Y m (u, v) = J K (2g + ^' ^ u s +-v s -, (21) 
y 4:71 (b — my (S + my. 

where mh is the projection of the angular momentum, —S < m < +S, u and v are the 
spinor components in spherical coordinates : 

u = cos(#/2)e^ /2 , v = sin(#/2)e^ /2 . (22) 

The radius R of the sphere is related to the number of vortices (or flux quanta in the 2DES 
language) that pierce it : 

R = (23) 

The one particle angular momentum S is such that 2S = N^. Due to sphere topology, 
the relation between the number of particles N and for a given fraction is linear with 
a non-zero shift. For each trial wavefunction for a given fraction, there is a unique shift, 
which is a characteristic of the quantum Hall state. In the case of parafermionic states, the 
relation between the magnetic flux and the number of particles is given by : 

N* = § N - 2 (24) 

This can be deduced from the expression of the on the sphere by applying a stereographic 
projection. Formally, we just have to drop the gaussian factor and make the substitution : 

(zi - Zj) — ► (uiVj - Ujvi) . (25) 

The two-body interaction is completely characterized by a set of 25 + 1 numbers {V m } 
called the pseudo-potentials^ 2 -. The integer m is the relative angular momentum between the 
two particles. For spinless bosons, only even-m potentials are relevant, s-wave scattering 
interaction corresponds to the case where all pseudo-potentials are equal to zero except 
Vo- Longer range interactions involve additional pseudopotentials, the next one for spinless 
bosons being V%. Thus, adding some Vi component allows to test the effect of longer range 
interactions. Comparison between the RR states or their quasihole excitations with the true 
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ground states is achieved by computing overlaps. For two states \^/) and |$), the overlap is 
defined as O = |($|\l/)| 2 . This definition can be extended to the case of subspaces of same 
dimension N and spanned by vector sets and {l^j)} : 

1 - 

o = ^£i<<m*;>i 2 - ( 26 ) 

Both exact ground states and RR/quasihole candidate ground states are evaluated us- 
ing exact diagonalizations of the associated Hamiltonian. Numerical diagonalizations are 
achieved using Lanczos-like algorithm or full diagonalization algorithm in a given L z sub- 
space. In the case of k + 1-body hardcore interaction, the matrix are being less and less 
sparse with increasing k value, requiring more memory and CPU time and making conver- 
gence harder to reach. Moreover due to the L z -only restriction on the Hilbert space, looking 
at the quasihole ground states require the evaluation of highly degenerate eigenstates. Thus 
we can reach lower system sizes compared to the ground state. 

V. GROUND STATE OVERLAPS 

We look at the overlap between the RR ground state and the exact ground state. Tables 
U] to IIVI display the overlaps for various fractions between the RR state and the two-body 
hardcore interaction hamiltonian ground state for different sizes. We also include overlaps 
with other ground states such as Coulomb interaction or n > 2-body hardcore interactions. 
In the particular case where N = 4z/, the overlap is equal to one. This is due to the dimension 
of the Hilbert subspace in the L = sector when S = 1 which is equal to one. 

Some of the results presented in table H] have already been published^. They show that 
the Pfaffian state is a good description of the physics at v = 1. As already noticed, longer 
range interactions tend to improve the overlap. 

For the v = 3/2,2,5/2 fractions, the situation is not so clear. Fewer values can be 
obtained and the overlap is non-monotonic with respect to the size of system, making dubious 
convergence to the thermodynamic limit. If we consider long-range interaction like Coulomb 
interaction, overlaps are improved, but we still get the same non-monotonic behavior. The 
same remarks are valid for the comparison with the n-body hardcore interaction (2 < n < 
k) : the overlaps are closer to unity as n tends towards k + 1. 
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To ascertain the role of longer-range interaction, we follow the method proposed in ref . (|12f ) 
for v = 3/2. We add a V2 contribution to the two-body hardcore interaction. Figure ^ shows 
the overlaps as a function of the ratio V2/V0 for the four filling factors v = 1, 3/2, 2 and 5/2. 
The conclusions we can draw are similar to ref. : long-range interactions help to stabilize 
the parafermionic ground state. Note that the drop of the overlap for large values of V2/V0 
is correlated to a similar effect in the gap value (see figure |2J) and is thus related to the loss 
of incompressibility. 



VI. QUASIHOLE EXCITATIONS 

If we believe that the parafermionic description is relevant for the fractions v = k/2, then 
quasihole excitations should also be present. Studying quasihole excitation on the sphere 
geometry is interesting on its own. Indeed, non-Abelian statistics is related to the quasihole 
ground state degeneracy. We can sort these states by their orbital quantum numbers L 
and L z . Evaluating the degeneracy of each sector is already a non trivial task. A formula 
was found for the Pfaffian case by Read and Rezayi^ 5 -. Gurarie and Rezayii^ have also an 
algorithm to compute the degeneracy in the v — 3/2 case. Finally, Ardonne has proposed 
an expression for the degeneracy valid for any v = k/2 value. We briefly describe how 
we extract the multiplet decomposition of the quasihole degenerate states from Ardonne's 
formula in an Appendix. Comparison of degeneracy values obtained using the CFT approach 
with the results of numerical exact diagonalizations is a way to validate the CFT approach. 
Numerical computations have been performed for the Pfaffian 15 at v — 1 and also^ for the 
v = 3/2 case. We give here additional values for these two fractions (tables IVl and IVI|) . We 
also compute degeneracies for v — 2 and v — 5/2 which haven't been published before (see 
tables IVTD and IV 111)) . The results we obtain are in agreement with Ardonne's formula. 

To test the validity of the quasihole hypothesis, we compute the overlap at a given fraction 
v = k/2 and for k quasiholes between the subspace spanned by the quasihole states of the 
/c + l-body hardcore Hamiltonian (jSJ) and the lowest energy states of the short-range problem 
at each L value. In each (L, L z ) sector, we thus consider the N^' q lowest energy eigenstates 
where N^ q is the degeneracy for q quasihole ground states at filling factor v = k/2 with 
angular momentum L (N^' q is L z independent) as candidates for the non-Abelian quasihole 
states. The corresponding overlap O k £ q is evaluated using definition Eq. (|26jl . In order to 
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easily characterize the agreement with the whole set of quasihole states for all L values, we 
introduce a total overlap defined as : 

ok , q = E L N k L q (2L + l)O k L « 
£ L i\£ 5 (2L + l) ' 

which is just another way to write the total overlap with respect to the subspace spanned 
by all quasihole states. 

Our results are given in tables ITXl IXl IXI| and IXHI for fractions v = 1,3/2,2 and 5/2 and 
q = k quasiholes. The success of the quasihole description is quite impressive at v — 1. 
However the agreement becomes increasingly worse with higher k values. Notice that for a 
given system with fixed k and N values, the smallest overlap is obtained for the largest L 
total momentum. Due to its high L z degeneracy, it plagues the total overlap. This certainly 
means that fewer quasihole excitation with non-Abelian statistics are present in the pure 
hard-core model than in the k + 1-body system. We can also add more quasiholes. Table 
IXIIII displays the results for q — 2k at v — 1. Considering the high degeneracy we are 
looking at (up to 336 for iV = 12), the overlaps are quite good especially if we do not take 
into account the ones associated to the largest total momentum. The effect of longer range 
interaction is similar to the ground state case : adding some Vi component tends to improve 
the overlap. In figureEl we have plotted the total overlap as a function of V2/V0 for fractions 
v = 1, 3/2, 2. There is now a maximum of the overlap for a moderate amount of longer-range 
interactions. 



VII. CONCLUSION 



Existence of quasiparticles with non-Abelian statistics is an exciting question of modern 
physics. The possible appearance of such quasiparticles in rotating ultracold boson gases 
is a strong motivation to experimentally reach the corresponding regime. At filling factor 
v = 1, we have shown that the pairing scheme of Moore and Read extends to the quasihole 
excitations. The degeneracy we observe is exactly that predicted by the CFT approach and 
the overlap of the subspaces spanned by the quasiholes shows that they are likely to be 
relevant at this fraction. Concerning the RR states for larger filling factors their overlaps 
with the ground state of the pure hard-core model are much less impressive and the size 
dependence is irregular, as was already observed in the gap values. We have found a set of 
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degenerate states with the quantum numbers predicted for quasiholes generated from the 
RR states by addition of flux quanta, they have also the features expected from the CFT 
approach. Again the overlaps (now for subspaces taken as a whole) are much less impressive. 
It is highly unlikely that the RR states are relevant for large k values. Adding long-range 
interactions like a second pseudopotential V2 certainly strengthen these RR states and their 
quasiholes states. So to construct in practice a RR state one may have to fine-tune the 
interaction potential between ultracold atoms. It is likely however that the Pfaffian at v — 1 
is the most conveniently implemented state for manipulation of non-Abelian statistics (but 
it does not support universal quantum computation^-). 
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IX. APPPENDIX 

Our purpose is to show how we can get the multiplet degeneracy for the quasihole states 
from Ardonne's formula. The formula gives access to the intrinsic degeneracy and is derived 
from the truncated characters of the algebra (or su(2)fc/u(l)). These truncated characters 
can be written as : 



k-l 

n 



f + (OU-i - en) .a), 

ft; 



(28) 



where the q-deformed binomial is defined as follow : 

nr=i (i - <? 



m 
P 



(29) 



n? =1 (i-^nr=7(w) 

and is equal to zero if p > m or m,p < 0. n is the number of quasiholes. It is linked 
to the number of added quantum fluxes A$ by the relation n = fcA$. Ik-i is the identity 
dimensional matrix and Ck-i — 2A fc : L 1 where Afc_i is the Cartan matrix of the su(k) algebra : 
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(Afe-i)ij = 28 i;j - 5\i-j\,i. (30) 

a = (oi, ...,ajfc_i) is a vector of — 1 non-negative integer such that Yli^i = F where 
F is a multiple of k. When we look at a system of N bosons, we are only interested in the 
values F = 0, k, 2k, N. F has to be understood as the number of unclustered bosons. 
In the simplest case k = 2, N — F corresponds to the number of parafermionic fields that 
appear when using OPE. Each g-polynomial in front of given x F monomial, is associated to 
the multiplet decomposition of the intrinsic degeneracy of a given F value. It will be of the 
form : 

El^ Et=-L <7° +L % (31) 

where a is a global shift of the q power and n\ is the number of multiplet of momentum 
Thus the multiplet decomposition of the intrinsic degeneracy can be directly read-out from 
the polynomial expression. 

For example, let evaluate (J28|) for k = 3, n = 6, For the N = 6 bosons case, we only need 
the partial development given in (f3~2J) up to x & : 

Y 6 (x; q,3) = 1 + x 3 (g 4 + q 3 + q 2 ) + x 6 q 6 + ... (32) 

Thus the multiplet decomposition for the intrinsic degeneracy is one singlet L = for F = 0, 
one multiplet L = 1 for F = 3 and one singlet L = for F — 6. The extrinsic part gives 
the following multiplet: one singlet L = for F = 6, one multiplet L = 3 for F = 3, one 
multiplet for L = 0,2,3,4,6 for F = 0. Using the standard momentum addition rules, we 
get the result display in table IVII 
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Figure 1: From left to right and top to bottom: overlap between the RR state and the ground 
state of longer range interaction Hamiltonian as a function of V2/V0 at filling factors v = 1, 3/2, 2 
and 5/2. 
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Figure 2: Gap of the longer range interaction Hamiltonian as a function of V2/V0 at filling factors 
v = 1 (a) and v = 3/2 (b). 



n ew Qi/r 

4 1.0 1.0 

6 0.9728 0.9728 

8 0.9669 0.9771 

10 0.9592 0.9659 

12 0.8844 0.9165 

14 0.8858 0.9213 

16 0.8833 0.9170 



Table I: Overlaps at v = 1 between the Pfafnan obtained as the ground state of the 3-body hardcore 
interaction hamiltonian and the ground state of the 2-body hardcore interaction or the Coulomb 
interaction (1/r) hamiltonian. 
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N O k=1 O k=2 1/r . 
6 1.0 1.0 1.0 

9 0.9642 0.9891 0.9642 

12 0.8647 0.9702 0.8904 

15 0.9189 0.9788 0.9307 

18 0.6774 0.9239 0.7226 



Table II: Overlaps at v = 3/2 between the RR state obtained as the ground state of the 4-body 
hardcore interaction hamiltonian and the ground state of the k + 1-body hardcore interaction or 
the Coulomb interaction (1/r) hamiltonian. 



Ok=l Ok=2 Ok=3 Ol/r 

8 1.0 1.0 1.0 1.0 

12 0.9636 0.9811 0.9949 0.9636 

16 0.7801 0.8919 0.9753 0.8037 

20 0.8822 0.9499 0.9874 0.8985 



Table III: Overlaps at v = 2 between the RR state obtained as the ground state of the 5-body 
hardcore interaction hamiltonian and the ground state of the k + 1-body hardcore interaction or 
the Coulomb interaction (1/r) hamiltonian. 



N Ok=l Ok=2 £>k=3 Ofc=4 Ol/r 

10 1.0 1.0 1.0 1.0 1.0 

15 0.9659 0.9789 0.9902 0.9975 0.9659 

20 0.7455 0.8291 0.9165 0.9798 0.7644 



Table IV: Overlaps at u = 5/2 between the RR state obtained as the ground state of the 6-body 
hardcore interaction hamiltonian and the ground state of the k + 1-body hardcore interaction or 
the Coulomb interaction {1/r) hamiltonian. 
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Figure 3: Upper part: total overlap for as a function of V2/V0 at filling factors v = 1 with two 
quasiholes (a, A<E> = 1) and four quasiholes (b, A$ = 2). Lower part: total overlap as a function 
of V 2 /Vq at filling factors v = 3/2 (c) and 1/ = 2 (d) for A$ = 1. 
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Table V: Number of multiplets of states at zero energy for the three-body Hamiltonian for v = 1 
for A<3? added quantum fluxes. # is the total number of degenerate states. 
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Table VI: Number of multiplets of states at zero energy for the four-body Hamiltonian for v = 3/2 
for A$ added quantum fluxes. # is the total number of degenerate states. L values for rows with 
a star, have to be understood as L — 1/2. 
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Table VII: Number of multiplets of states at zero energy for the five-body Hamiltonian for v = 2 
for A$ added quantum fluxes. # is the total number of degenerate states. 
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Table VIII: Number of multiplets of states at zero energy for the six-body Hamiltonian for v = 5/2 
for A<I> added quantum fluxes. # is the total number of degenerate states. L values for rows with 
a star, have to be understood as L — 1/2. 
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N 


2,2 


Co' 2 


of 


of 


of 


of 


of 

•J 


of of 


4 


1.0 


1.0 




1.0 










6 

8 


0.8579 
0.8760 


0.9857 


0.9994 


0.8884 


0.7972 


0.8569 






10 


0.8661 




0.9111 




0.8334 




0.8287 




12 


0.6800 


0.9651 




0.8007 




0.8392 




0.5015 


14 


0.6825 




0.7685 




0.7059 




0.7417 


0.611 



Table IX: Overlap of the lowest energy exact wavefunctions and the Pfaffian two quasihole wave- 
functions at v = 1. 



N 



O 



3,3 



o\ 



3,3 



o\ 



3,3 



of 



O 



3,3 



Of O 



3,3 
5 



of of 



12 
15* 



1.0 
0.7933 
0.6289 
0.4440 



0.8221 



1.0 

0.9561 

0.7965 



0.7065 
0.7317 
0.5754 



1.0 

0.8579 
0.5856 



0.7803 

0.5735 0.4897 

0.4440 0.7572 0.0009 



Table X: Overlap of the lowest energy exact wavefunctions and the k = 3 three quasihole wave- 
functions at v = 3/2. L values for rows with a star, have to be understood as L — 1/2. 
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0.0957 



Table XI: Overlap of the lowest energy exact wavefunctions and the k = 4 four quasihole wave- 
functions at v = 2. 
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e> 5 ' 5 


^0 
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O 5 / 
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0.8110 




0.2827 



Table XII: Overlap of the lowest energy exact wavefunctions and the k = 5 five quasihole wave- 
functions at v = 5/2. L values for rows with a star, have to be understood as L — 1/2. 



N O 2 ' 4 O 2 A 0\ A 0\ 



o 



2,4 



o 2 / ol A 



o 2 > 4 o 2 / ol A 



o 



2.1 



ro 2,4 , 2,4 , 2,4 
^10 ^11 u 12 



4 0.9807 1.0 
6 0.8541 0.9998 
8 0.8458 0.9036 
10 0.6692 0.8991 
12 0.5788 0.7964 



1.0 0.9572 
0.9741 0.9764 0.9369 0.5590 
0.9407 0.9197 0.8155 0.7370 0.8620 0.8184 
0.8269 0.7729 0.8281 0.7291 0.8090 0.4411 0.5091 0.2898 
0.7970 0.8253 0.7832 0.7756 0.5882 0.5654 0.5572 0.4255 0.4489 



0.0692 



Table XIII: Overlap of the lowest energy exact wavefunctions and the Pfaffian four quasihole 
wavefunctions at v = 1. 
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